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Abstract 

1^ . 

The necessary and sufficient condition for a quantity to be mea- 
surable by the HBT method is given and discussed. 



> 

(N 1 Introduction 

(N 
O 

^ I The HBT method gives valuable information about the interaction regions 

Q ' in multiparticle production processes, but it has important limitations. For 

^ ■ instance, it is well known, cf. e.g. the review [1] and references quoted 

there, that the total interaction region cannot be reconstructed. The best 
one can hope for, is the reconstruction of the separate homogeneity regions. 
^ ' Homogeneity region K means the region of space, where the particles with 

' momentum K are produced. The geometry of homogeneity region K can 

be described by specifying the profile function p(X|K), i.e. the probability 
density of the production points for the particles with momentum K. Instead 
of the profile functions one could use their moments, or central moments and 
(x)(K), or cumulants. 

In the present paper we address the question: how much can be learned 
about the profile functions using the HBT method. The use of cumulants 
will be particularly convenient. 
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2 Simplifying assumptions 



In order to make the problem tractable and/or to avoid unessential compli- 
cations we make three simplifying assumptions. 

• The experimental data is perfect. This means that for any n the n- 
particle momentum distribution can be measured with arbitrary preci- 
sion and fully corrected for final state interactions, resonance produc- 
tions etc. 

• The profile function can be related to the emission function ^(X, t, K) 
by the formula 

p(X|K) = JdtS{X,t,K), (1) 

where the emission function can be interpreted both as a classical, time 
dependent phase space distribution, then the relation is obvious, and 
as a generalized Wigner function. This assumption about the interpre- 
tation of the emission function is made in most analyses of the HBT 
method. Rough estimates suggest [2] that the error introduced by it 
is negligible for heavy ion collisions and probably also acceptable for 
e'^e~ annihilations. The interpretation as generalized Wigner function 
is necessary to connect the emission function, and consequently the pro- 
files, to the measured single particle density matrix in the momentum 
representation. 

• The momentum distribution for n identical bosons can be expressed 
in terms of an effective, n-independent, single particle density matrix 
p(p; p') as follows [3] 

n 

p(Pi, . . . , pj = c„ ^ n p(p.; pp.)> (2) 

p j=l 

where C„ are normalizing constants and the summation is over all the 
permutations of the momenta p^, . . . , p,„. Thus, e.g. the single particle 
momentum distribution is p(p; p) and for the two-particle distribution 
we get the well-known formula 

^(Pi,P2) ~ p(Pi;Pi)p(P2;P2) + |p(Pi;P2)P- (3) 
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In the absence of final state interactions, resonances etc. tliis is a 
standard assumption. 

It can be made plausible as follows. P{p^, . . . ,p^) is the diagonal 
element of the n-particle density matrix in momentum representation. 
For independent particles 

n 

PiPl,---,Pn)^I[piPpPj)- (4) 

This corresponds to the well-known result of the theory of probabilities 
that for independent events the probability of a set of events is equal to 
the product of the probabilities of the single events. Bose-Einstein sym- 
metry implies that the n-particle density matrix p(pi, . . . , p„; , . . . , p^) 
must be invariant with respect to permutations of p^, . . . , p„ and with 
respect to permutations of p'^, . . . , p^. Symmetrizing (4) with respect 
to the second arguments one obtains (2). This symmetrization makes 
the result symmetric also with respect to permutations of the first argu- 
ments, but it spoils the normalization for n > 1. Thus, the normalizing 
factors Cn have to be introduced. It is seen that assumption (2) can 
be interpreted as the simplest symmetrized version of (4). 

3 Ambiguities 

It is easily seen that the right hand side of formula (2) is invariant with 
respect to the substitution 

p(p;p')^p'(p;p') = p(p';p)- (5) 

What is more interesting, it is also invariant with respect to the substitution 
[4] 

p(p;p') ^ p'(p;p') = e^^(P)p(p;p')e-^(P'), (6) 

where /(p) is an arbitrary, real-valued function of p. Transformations (5) 
and (6) generate the full invariance group of expression (2) [5]. Thus, the 
problem of finding the uncertainties of the deduced profile functions reduces 
to the problem of finding the implications of the ambiguities (5) and (6). 
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4 Cumulants 



Let us introduce the familiar notation 



K = - 




q = Pi - P2- 



(7) 



Then, using our second simphfying assumption, i.e. interpreting the emission 
function as the generahzed Wigner function: 



Using the terminology of the probability calculus this means that ^^-^Qj at 

given K is the characteristic function of the probability distribution (profile 
function) p(X|K). 

The cumulants of profile functions will be denoted ]C{rx,ry,rz), where 
i^x,i^y, fz are non- negative integers and the argument K is not written explic- 
itly. The number r = r,j. + Vy + will be called the order of the cumulant. 
The cumulant is even (odd) when r is even (odd). The cumulants are related 
to the characteristic function by the standard formula 



Expanding the exponent in the integrand in formula (9) one finds that 



with similar formulae for the other two order one cumulants. These are the 
only cumulants which cannot be expressed in terms of the central moments 
of the profile function and thus give information about the absolute posi- 
tions of the homogeneity regions. The order two and three cumulants form 
respectively the covariance matrix and the skewness matrix. The fourth or- 
der cumulants are more complicated, because they depend on both second 
order and fourth order central moments. This is easily understood because, 
as easily checked, for a Gaussian probability distribution all the cumulants 
of order higher than two vanish. 







/C(1,0,0) = {x) 



(10) 
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5 Result 



Let us consider the effect of the ambiguities (5) and (6) on the deduced profile 
functions. Ambiguity (5) means that the data cannot distinguish between 
the effective single particle density matrices p(K, q) and p(K, — q). This 
corresponds to the space inversion of the distribution p(X|K). Since space 
inversion changes neither the size nor the shape of the interaction region, this 
ambiguity is not very interesting. 

Ambiguity (6) means that the data do not distinguish between p(K, q) 
and p'(K,q) where: 



Expanding the logarithms on both sides according to (9) it is immediately 
seen that the ambiguity affects only the odd cumulants. Explicitly, the rela- 
tion between the cumulants is: 



2 ) d^y d"' 

For all real valued functions /(K) this formula gives the full ambiguity for 
the cumulants and consequently also for the profile functions. Thus, it solves 
our problem. 



6 Discussion 



The even cumulants {r — rx-\-ry-\- even) can be unambiguously measured 
[6]. In particular: 

• The HBT radii, which are given by the r = 2 cumulants, can be unam- 
biguously determined [7]. 

• The kurtosis matrix, which can be expressed by the r = 4 and r = 2 
cumulants, can be unambiguously measured. The same holds for all 
the even central moments of the profile functions p(X|K). 
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• The distribution of Xi — X2 for pairs of points in the homogeneity 
region, which is important for the imaging method [8], [9], [10], [11], 
can be unambiguously measured, because it can be expressed in terms 
of even cumulants only [6]. 

The odd cumulants, which are as important for the determination of the 
interaction region as the even ones, are not measurable. 

In particular, the relative positions of the centers of the homogeneity 
region are almost unconstrained. Since 

(x)(K) = V/(K) (13) 

and /(K) is arbitrary, it is enough to assume that (for sufficiently well- 
behaved functions /(K)) the rotation of the field (x)(K) vanishes. This 
explains, for instance, the behavior of the exploding source model described 
in [12]. There, the known input-size of the interaction region is different 
from the size which would be obtained by analyzing the momentum distri- 
bution under the assumptions currently used by experimentalists. In our 
terminology, the function /(K) assumed in the input is different from the 
corresponding function tacitly assumed in the analysis based on the momen- 
tum distributions. Consequently, the homogeneity regions get shifted with 
respect to each other and the overall interaction region changes. In general, 
given a momentum distributions for 1,2,3,... identical bosons which cor- 
respond to a reasonable interaction region, it is easy to construct a model 
where exactly the same momentum distributions correspond to an interaction 
region the size of a football. 

On the other hand, the same function /(K) describes the ambiguities for 
all the odd cumulants. Thus, shifting the homogeneity regions one usually 
has to deform them. Consider the following example [4]. Assume that the 
emission function is proportional to 5{t). Then the integration over time 
yields the ordinary Wigner function. We start with the model where 

1 K 1 D2r»2 

p(Pi,P2) = ^e-iK--^«P. (14) 

and are positive constants. The corresponding profiles are all Gaus- 
sian, centered at x = 0. Then we deform the density matrix according to (6) 
with 
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,3 ' 



(15) 



where Qj arc positive constants. This model can be solved analytically and 
leads to non-trivial observations. For aj — > oo there is no correction to the 
original Gaussians. When the parameters decrease a number of things 
happens. The centers of the homogeneity regions shift so that 



This increases the overall interaction region. The main peak in p{xj\K) be- 
comes steeper. To the left from it secondary maxima develop. Between these 
maxima there are regions where the Wigner function is negative. Negative 
values of the Wigner function are inconsistent with its interpretation as a 
classical phase space density. This is a well-known problem (see e.g. the 
reviews [13], [14]) which is usually solved by smearing the Wigner function. 
A suitably smeared Wigner function can be interpreted as a phase space den- 
sity, but it is not the original Wigner function any more. These are problems 
beyond the scope of the present report. 

The odd cumulants, though separately unmeasurable, are strongly corre- 
lated with each other. If the distribution (x)(K) is known, one can use equa- 
tion (13) to find, up to an irrelevant constant, function /(K) and then all 
the odd cumulants are unambiguously determined. Since the even cumulants 
arc measurable, this implies that the profiles p(X|K) are also unambiguously 
determined. Models based on macroscopic physics, like the hydrodynamic 
models, are more likely to give correctly the distribution (x)(K) than the 
finer details of the shapes of the interaction regions. The fact that getting 
(x)(K) is enough to determine unambiguously from experiment the profile 
functions may thus be of practical importance. 
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